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INTRODUCTION 


The  radiation  produced  by  gamma  rays  incident  on  ordinary 
dielectric  materials  such  as  glass  was  first  discovered  by 
Cerenkov^  in  1934  and  was  described  in  terms  of  a  charged 
particle  (electron)  moving  faster  than  light  in  the  medium  by 
Frank  and  Tamm^  in  1937.  A  summary  of  work  to  195Q  is  contained 
in  the  treatise  by  Jelly 3.  An  important  application  is  the 
Cerenkov  particle  detector  which  is  fami?  Lar  in  any  particle 
physics  laboratory,  and  an  early  and  crucial  application  occurred 
in  the  discovery  of  the  antiproton4. 

Because  the  distribution  of  intensity  of  Cerenkov  radiation 
is  proportional  to  the  frequency,  the  radiation  at  microwave 
frequencies  would  be  low  unless  beams  are  intense  and  bunched  so 
that  coherent  radiation  by  many  electrons  contributes.  Danos5 
in  1955  calculated  radiation  produced  by  a  planar  beam  passing 
above  a  dielectric  interface  and  a  hollow  cylindrical  beam  passing 
through  a  hole  in  a  dielectric.  Experimental  and  theoretical 
investigations  at  microwave  frequencies  ware  reviewed  by 
Lashinsky6  in  1961. 

This  investigation  was  motivated  by  a  recent  renewed  interest 
which  has  included  the  study  of  stimulated  Cerenkov  radiation,  in 
which  the  electron  may  be  in  a  medium  consisting  of  a  gas7  or  a 
hollow  dielectric  resonator8'9.  Recent  developments  of 
electron  accelerators  for  applications  such  as  free  electron 
lasers  (PEL)  have  aimed  toward  high  peak  currents  in  bunches  in 
contrast  to  nuclear  and  particle  physics  applications,  where  low 
peak  but  high  average  currents  are  desirable  to  avoid  saturating 

1 


detectors.  The  high  peak  currents  in  the  new  accelerators  should 


j  .a  id 


enhanced  Cerenkov  radiation,  as  is  calculated  in  this 


paper . 


CALCULATION  OF  THE  POYNTING  VECTOR 
In  the  following  derivation,  we  consider  the  Cerenkov 
radiation  produced  in  a  dispersionless  medium  such  as  gases  or 
other  dielectrics,  by  a  series  of  pulses  of  electrons  such  as  ar 
produced  by  a  traveling  wave  electron  accelerator  (Linac).  The 
pulses  or  bunches  are  periodic,  the  total  emission  region  is 
finite  and  the  bunches  have  a  finite  size. 

In  determining  the  radiated  power,  the  procedure  is  to 
calculate  the  Poynting  vector  from  fields  which  are  in  turn 
obtained  from  solutions  of  the  wave  equations  for  the  potentials 
Since  the  current  and  charge  densities  entering  into  the  wave 
equations  are  expressed  in  fourier  form  the  resulting  fields  and 
radiated  power  also  have  fourier  components.  In  the  derivation, 
is  the  coordinate  at  which  the  fields  will  be  calculated,  r/  is 
the  coordinate  of  an  element  of  the  charge  which  produces  the 
fields  and  n  is  a  unit  vector  in  the  direction  of  iT.  We  assume 
that.  E(r,t.)  and  B(r,t.,)  have  been  expanded  in  a  fourier  series, 
appropriate  for  the  case  where  the  source  current  is  periodic - 
Then  we  have 

00 

E  Cr,t)  =  ^  e  (x,a>) 

co=—  (1) 

and  a  corresponding  expansion  for  B,  where  a>  is  a  discrete 
variable  and  E  and  B  are  fourier  series  coefficients.  The 
poynting  vector  S  is  given  by 


and  it.  is  easy  to  show  that  the  average  of  S  in  a  direction  given 

by  a  normal  vector  n  is 
T  <» 


1 

T 


O 


A  ^ 

n-S 
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(3) 


where  T  is  an  integer  multiple  of  the  period  of  the  periodic 
current. 

Letting  c  =  (ye) ”1/2  be  the  velocity  of  light  in  the 
medium,  the  wave  equations  for  A,4>  and  their  solutions  are, 

(v  --2— 2]  A(r,t)  =  yj(r,t) 

»  c  dt  /  (4) 

(  v  —  ~2^ — 2 J  ~  i/2  p(r,t) 


A(r, t) 


u  t~t ’)  J(r \t ’)  d3r  'dt ' 

1 JJJJ  D  (r-r  ’  ,  t-t' )  p  (r'  ,t’)d3r'dt' 


(5) 


where  the  Green's  function  D  is  given  by 


D(r,t)  =  _ 1_  6 (t-r/c) 

4irr 


(6) 


The  vector  potential  A(r,  t.)  also  can  be  developed  in  a  fourier 

series  expansion  of  a  form  similar  to  (1)  with  an  expression  for 

the  fourier  series  coefficients  given  by 

-  1  f  dtA(r,t)eiut 
T  J 

o 


4Tr  I  r-r1 1 


/c 


(7) 


Now  if  we  assume  that,  the  observer  is  far  from  the  source  so 
that  fr j>> Jr*  J  for  regions  where  the  integrand  in  (7)  is  important 
we  can  let  |r  -  r '  |  =  r  -  n  •  r '  in  the  exponential  and  Jr*  -  r*  |  = 
in  the  |r  -  r'  |“1  factor  in  (7),  obtaining  (where  n  =  r/r) 

A{r,w)  -  u  e1"170  ,rr/'d3r'7(?',u)e-i(“/c,fi-?' 

4irr  JJJ 


The  fourier  series  coefficients  of  the  fields  are  obtained 

from  those  for  the  vector  potential  (8)  through  the  usual 

relations  B  =  "7  x  A  and  £  =  -  7<J>  -3A  .  Under  the  conditions 

3 1 

leading  to  (8)  the  field  fourier  coefficients  are^; 


^  A  -k  ^ 

B,(r,0))  =  ico  n  x  A(r,<a) 
c 

E(r,co)  =  -c  n  x  B  (r,cu) 


(9) 

(10) 


The  poynting  vector  can  now  be  found  by  using  (9)  and  (10) 

in  expansions  like  (1)  and  then  substituting  in  (2).  However  it 

is  more  convenient  to  deal  with  the  frequency  components  of  the 

radiated  power  by  substituting  (9)  and  (10)  into  the  expression  of 

the  average  radiated  power  (3) . 

T  » 

f  n*Sdt  =  1  o2  |  n  x  A(r,o>)  j  2 
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FOURIER  COMPONENTS  OF  THE  CURRENT 
The  expression  (7)  for  the  fourier  components  of  the  vector 
potential  contains  the  fourier  components  of  the  current  density. 
Consequently  it  is  necessary  to  examine  the  form  of  the  current 
and  its  fourier  development.  Assume  the  current  is  in  the  z 
direction  and  periodic.  If  the  electrons  move  with  velocity  v, 
and  we  ignore  for  the  moment  the  x  and  y  variables,  the  charge  or 
current  functions  should  have  the  general  form 


ik  Z  _  -itot 


(z,t)  =  £  e  z  £ 


e  "z  's  e  ~  f(k  ,co) 

z 


01 


Under  the  condition  of  rigid  motion, 


ffz,t)  =  f0(z-vt) 


(12) 


(13) 


it  is  easy  to  show  that 


k  vi,<V 

z 


(14) 


where 


f  (kj 

-v  a  Z 


1 

z 


e  XKzz  f0(z)dz 


(15) 


Thus  the  restrictions  of  equation  (13)  reduce  the  two 


dimensional  fourier  series  of  eq.  (12)  to  essentially  a  one 
dimensional  series  (14). 


With  (14)  in  mind,  the  current  density  associated  .it h  the 
electron  beam  from  a  linear  accelerator  should  be  periodic  in  bet 
z,  t,  with  a  fourier  series  expansion,  but  the  x  and  y  dependence 
should  be  represanted  by  a  fourier  integral  form: 


+00 


Jz(r,t)  =  vp(r,t: 


=  — — _  f  dk  fdk  Y”  e  ^ 

(2,)2  }  x  J  y  L  _ 


(k  •  r  -  (A  t ) , 


-00  -oo 


kz=- 


16) 


where  the  fourier  components  of  the  charge  density  are 

to  00 

=  / dx  /  dy  k 

—  cc  -  CO 

pQ  C r )  is  p(r\t)  evaluated  at  t  =  o  and  J  is  assumed  to  be  in 
the  z  direction.  Note  in  eq.  (16)  that  kz  and  w  are  both 
discrete  and  from  (14),  a)  =  kzv. 


U 

I 


dze 


-ik  ■ 


rp0 (r) 


(17) 


VECTOR  POTENTIAL 


The  results  of  the  previous  section  can  be  applied  to  the 
evaluation  of  the  vector  potential  and  in  turn  to  the  fields. 

Let  the  infinite  periodic  pulse  train  be  made  finite, 
extending  from  z  =  -2 ‘  to  z  =  +Z  1  and  let  9  be  the  angle  between  n 
and  A.  Then  the  cross  product  in  (11)  can  be  written 


■?/**■  »  i  .  a  y  icur/c 

n  x  A(r,co)|  =  sm9  e 


®  CD  Z 1 

/  dx'  /  dy '/ 

sco  xoo  —7,  1 

^  2 


TSFJ 


dz  1 


-in-r  'o)/c 


dk 
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dkY 


w 
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k2=_» 


vD-°<k)5k  v,»e 


ik  •  r ' 


(18) 


But 

r  7-5'  a 

f  dx'  I  dy'  /  do'  elr'-,k  -  ™/o) 

-i.®  -Z  ' 

=  ( 2tt)  2  5(kx  -  nx0)/c)  <S  (ky  -  nya)/c)I(Z')  (19) 


where 


KZ 


J-7.' 


dz'  ei<kz  -  nz“/c)*'  .  2  sinG2. 

G 


and  G  =  kz  -  nz  oj/c  *  oj/v  -  nz  oj/c 
And  thus  the  cross  product  term  is 


(20) 


I  n  x  A(r,o>)  |  »  sin04i7  eicor/,c  v£o(nx'o/c,  n  oj/c,  co/v)  I  (Z  '  ) 


Note  that  u  is  a  discrere  variable  but  from  19,  the 


continuous  variables  kx  and  ky  become  evaluated  at  discrete 
points . 

Returning  to  (17),  a  more  symmetric  form  may  be  obtained  by 
assuming  that.  pQ(r)/  which  is  periodic  in  z  with  period  Z,  is 
actually  zero  between  the  pulses.  Denoting  by  pQ ' (  r )  the  charge 
density  of  a  single  pulse,  which  is  zero  for  z  <  o  and  z  >  Z  the 
integral  on  z  can  be  written 

00 

,  -ik  z  f  ,  -ik  z  ' 

dz  e  z  p0(r)  =  dz  e  z  pa  (r) 

J  (22) 

0  O  — 00 


J 


dz 


e“ikzzp6(r)  = 


■/ 


Then  (17)  the  fourier  coefficient,  of  the  charge  density,  becomes 


(23) 


where£0'(k)  is  the  three  dimensional  fourier  transform  of  the 
single  pulse  decribed  by  p  0 ' (r) . 

Substituting  these  expressions  into  (21)  gives  a  final  simple 
result  for  the  cross  product,  form: 


where 


I(Z'  )  sin  GZ1 

G  =  u/v  -  n  oj/c  (25) 

z 

k  =  (n  ai/c,n  (ii/c,ai/v) 
x  y 


The  components  of  the  Cerenkov  E  and  B  fields  may  now  be  found  by 
substituting  (24)  in  (9)  and  (10). 
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RADIATED  POWER 


The  frequency  components  of  the  average  radiated  power  are 
obtained  by  substituting  (24)  into  (11).  The  negative  frequency 
terms  equal  the  corresponding  positive  frequency  terms,  yielding 
factor  of  2  when  the  summation  range  is  changed.  Multiplying  by 
r2  converts  to  average  power  per  unit  solid  angle,  dP/dfl, 
yielding 


(26) 


G 


a 


where  w(o)*n)  is  defined  to  be 
2 

W(<o,n)  *  ^  ,  ff-  sin20(v2/Z2:  |o0*  (k)  |2  I2(Z') 
(4tt)  2  C 


(27) 


W(o>,n)  is  the  power  per  unit  solid  angle  radiated  at  t.ho  frequency 
co,  which  is  a  harmonic  of  the  basic  angular  frequency  co0  of  the 
periodic  pulse  train. 

To  find  ,  the  total  power  radiated  at  the  frequency  at , 

W  is  multiplied  by  d!2  and  integrated  over  solid  angle.  dote  that 
nz  *  co 8  9,  and  as  0  varies,  G  changes  according  to  (25), 


with  dG  «  -  (<o/c)  dn  so  that 

z 


dft  =  d$  ( c/ co )  dG 


Noting  that  the  integral  over  $  yields  2 * ,  we  find  the  result  for 
the  total  radiated  power  at  the  frequency  u>  for  all  angles 


fe|2  4  /  sin^9|gl(k)|2  I2(2')  w  dG 


lj»*  A  **  ’»  »  *•  %  »  '  »  ’  *  •  <*  •  *  •  .  *  »  *  v*  »•«»,* 


CERENKOV  ANGLE 

The  remaining  integral  over  G  may  now  be  examined.  The 
sin^  9  and  pQ  factors  may  often  be  slowly  varying  compared  to 
the  l2(z*)  factor,  the  latter  being  shown  in  Fig.  1.  For  large 
Z',  the  peak  in  I ^ ( Z  * )  becomes  narrow,  and  if  the  integrand  may 
be  neglected  outside  the  physical  range  G' <G<G" , 


r2  I 


I  (Z')dG  =  f  4  ( Z  ) 


i.2  I sinGZ  \ 2 


‘dG  =  4tt2  ' 


Then,  evaluating  the  sin  9  factor  and  p ' (k)  at  the 
point  corresponding  to  G  *  0,  (which  is  cos  9  -  nz  =  c/v)  shows 
that  0  at  the  peak  of  I(Z')  is  the  usual  Cerenkov  angle  9C.  We 


thus  obtain  for  large  Z ' 


P,t  -  xr  wv2sin20  |  g0  (k)  |  2  4irZ  '/Z 


Now  let.  2Z'/Z  *  ratio  of  the  interaction  length  to  pulse  spacing  = 
N,  the  number  of  pulses.  Also  S  *  v2  i/  at,  or  2tt  /  %  =  0)o/v  30 
that,  (in  the  large  Z'  limit), 

P0)  ’  47  ^•vsin20c|p;(ic)  |2N. 

(32) 

To  compare  with  usual  formulations,  (32)  is  divided  by  Nv  to 
obtain  the  energy  loss  per  unit  path  length  per  pulses 


If  *  JF  single;  (if)!2 


If  the  pulse  is  in  fact  a  point  charge,  the  fourier  transform 
Oq  (k)  reduces  to  q,  the  total  charge  per  pulse  and  (33)  is 
very  similar  to  the  usual  Cerenkov  energy  loss  formula,  where  for 
a  single  charge  q,  the  radiation  is  continuous  and  the  factor 
oj  in  (33)  is  replaced  by  uido.  In  the  present  case  the  pulse 

train  is  periodic  at  angular  frequency  oi0  and  the  radiation  is 
emitted  at  the  harmonic  frequenies  denoted  by  co  . 


DISCUSSION  OF  RESULTS 


Equation  (29)  and  the  approximate  evaluation  expressed  as 
(32)  form  the  main  results.  Some  consequences  will  now  be  noted. 

a.  EFFECT  OF  PULSE  SIZE.  The  spatial  distribution  of  the 
charge  in  the  pulse  appears  in  £q(£),  is  the  fourier 

transform  of  the  charge  distribution.  The  peak  of  in 

figure  1  occurs  at  G  -  0  or  nz  =  c/v.  Thus  at  the  peak,  a/v  =  n2  , 

so  that  k,  the  argument  of  p*  (Jc) ,  is  evaluated  at  (from  25) 

**  0 

k  *  nm/c 

(34) 


We  may  also  define  a  charge  form  factor  F(k) 

P*  (K)  -  qF(k) 

—  J 


The  form  factor  F(k)  is  identically  one  for  a  point  charge, 
and  for  a  finite  distribution  F(1c)  *  1  for  k  *o. 

Furthermore  F(k)  must  fall  off  as  a  function  of  Tc  near  the 
origin  if  all  the  charge  has  the  same  sign.  If  the  pulse  were 
spherically  symmetric,  F(k)  would  behave  as  elastic  electron 
scattering  form  factors  defined  for  nuclear  charge 
distributions11.  In  that  case,  the  mean  square  radius  <r^> 
of  the  charge  distribution  is  given  by  the  behavior  of  F(k)  near 
the  origin. 


F(k)  •*  1  -  <r^>  k^/6  (spherical  pulse)  (36) 


b.  SMEARING  OF  THE  CERENKOV  ANGLE.  For  a  finite  region  over 


which  emission  is  allowed,  namely  if  2Z '  is  finite,  the  function 
I2(Z'),  appearing  in  the  integral  in  (29),  will  have  a  finite 
width.  Since  the  peak  height  is  4Z'2  and  the  area  is  4^Z', 

(30),  we  can  assign  an  effective  width  2T  =  area/  height  *  tt/ z ’  * 
or 

T  =  ir/2Z*  (37) 

Thus  the  radiation  is  emitted  mainly  near  G  =  o  (which 
corresponds  to  9C)  but  in  a  range  Ag  »  +T  .  But  from  (25), 

0)  CO  -  x 

AG  -  —  An-  *  —  a  (cos 9)  so  that  there  is  a  range  in  cose  over 
c  c 

which  emission  occurs: 


A  (cosQ  )  *  £  ^r 


(38) 


Note  that  the  finite  angular  width  of  the  Cerenkov  cone  angle 
in  (38)  has  the  factor  l/o» ,  indicating  that  the  higher  harmonics 
are  emitted  in  a  sharper  cone. 

c.  BEHAVIOR  AT  HIGH  FREQUENCIES  RELATED  TO  PULSE  PARAMETERS. 
To  be  specific  let  the  charge  distribution  for  a  single  pulse  be 
given  by  gaussian  functions 


p0(r)  ■  A  exp(-x2/a2  -  y2/a2  -  z2/b2) 


(39) 


Then  F(k)  may  be  found  : 


F(k)  ■  exp(-k, 


,V/4 


-k 


V/4 


-k 


V/4, 


Beam  pulse  parameters  could  then  be  determined  by  measuring  the 
Cerenkov  radiation.  For  example,  fast  electrons  from  an 
accelerator  in  air  will  emit  with  a  9C  of  several  degrees 
in  which  case  kx  and  k„  in  (40)  can  be  neglected,  giving 

F(k)  »  exp(-kz2b2/4)  =  exp  [-u2b2/  (4v2)  ]  (41) 

The  expected  behavior  of  P  as  a  function  of  co  is  shown 
qualitatively  in  Fig.  2  as  a  linear  rise  at  low  frequencies 
followed  by  a  fall  off  at  higher  frequencies,  the  peak  occurring 
at 

<d  =  v/'b 
m 

(42) 

Furthermore,  a  different,  behavior  would  be  expected  at  very 
high  frequencies.  The  formulation  from  the  beginning  represents 
coherent,  radiation  from  all  charges,  not  only  from  one  pulse,  but 
coherence  from  pulse  to  pulse.  F(k)  then  describes  interference 
of  radiation  emitted  from  different  parts  of  the  pulse,  but.  note 
that  expressions  (29)  and  (32)  will  still  be  proportional  to 
q2  *  n^e^  where  n  is  the  number  of  electrons  in  a  pulse. 

Thus  the  n^  dependence  of  P  indicates  coherence.  But.  above 
some  high  frequency  such  that  w^/c  *  2ir /l,  where  l  is  the 
mean  spacing  of  electrons  in  the  cloud,  the  radit.ion  should  switch 
over  to  incoherent  radiation  from  each  charge  and  P^  should  be 
proportional  to  n.  The  incoherent  radiation  should  then  rise 
again  as  a  function  of  u> . 


CONCLUDING  REMARKS 


The  general  results  presented  here  describe  the  Cerenkov 
radiation  produced  by  fast  electrons  produced  by  a  linear 
accelerator.  For  an  S  band  Linac  operating  at  about  3Ghz  (10  cm 
radiation) ,  the  electron  bunches  are  separated  by  10  cm  and  would  be 
about  1  cm  long  at  1%  energy  resolution.  Microwave  Cerenkov 
radiation  is  expected  and  has  been  seen  in  measurements  at  the  Uaval 
Postgraduate  School  Linac. 

Two  types  of  measurements  were  made.  In  measurements  of  Series 
A,  an  X-band  antenna  mounted  near  the  beam  path,  oriented  to 
intercept  the  Cerenkov  cone,  was  connected  to  a  spectrum  analyser. 
Harmonics  3  through  7  of  the  2.85  GHz  bunch  frequency  were  seen  but 
power  levels  could  not  be  measured  quantitatively.  Harmonics  1  and  2 
were  below  the  wave  guide  cut  off.  In  the  series  B  measurements,  the 
electron  beam  emerged  from  the  end  window  of  the  accelerator,  a.id 
passed  through  a  flat  metal  sheet  90  cm  downstream  oriented  at  an 
angle  <t  from  the  normal  to  the  beam.  The  metal  sheet  defined  a 
finite  length  of  gas  radiator,  and  reflected  the  Cerenkov  cone  of 
radiation  toward  the  accelerator  but  rotated  by  an  angle  2$  from  the 
beam  line.  A  microwave  X-band  antenna  and  crystal  detector  with 
response  from  7  to  above  12  GHz  could  be  moved  across  the  (reflected) 
Cerenkov  cone  as  a  probe. 

As  mentioned  earlier,  the  series  A  measurements  showed  the 
radiation  is  produced  at  the  bunch  repetition  rate  and  its  harmonics. 
Series  B  measurements  performed  with  several  antennas  always 
indicated  a  broadened  Cerenkov  cone  with  strong  radiation  occuring  at 
angles  up  to  10*,  well  beyond  the  predicted  Cerenkov  angle  of  1.3*. 


Since  a  broad  band  detector  was  used  it  was  impossible  to  verify 
the  prediction  (see  eq .  38)  that  the  broadening  of  the  cone  should 
depend  on  the  harmonic  number.  However,  it  should  be  r.oted  that 
incoherent  radiation  by  a  beam  of  lu  A.  at 9  c  =  1.3°  for  a  1 
meter  path  in  air  would  be  about  10“14  watts  at  microwave 
frequencies  so  that  observation  of  any  signal  by.  either  method  A 
or  B  clearly  demonstrated  coherent  radiation  by  the  electron 
bunches . 

Many  of  the  concepts  were  clearly  noted  by  Jelly  in  his 
treatise  (Jelly^,  Section  3.4  especially).  The  form  factor  was 
noted  but  a  general  expression  was  not  given.  In  fact,  the  form 
factor  quoted  by  Jelly  represents  the  special  case  of  a  uniform 
line  charge  of  length  L'  with  a  projected  length  L-L‘ccs9c  in 
the  direction  of  the  radiation.  The  coherence  of  the  radiation 
from  the  bunch  was  noted  but  no  broadening  of  the  cone  nor 
harmonic  structure  were  developed. 

Casey,  Yeh  and  Kaprielian-^  considered  an  apparently 
related  problem  in  Cerenkov  radiation,  in  which  a  single  electron 
passes  through  a  dielectric  medium,  where  a  spatially  periodic 
term  is  added  to  the  dielectric  constant.  The  result  is  radiation 
occurring  even  for  electrons  which  do  not  exceed  the  velocity  of 
light  in  the  medium,  and  at  angles  other  than  the  Cerenkov  cone 
angle.  The  non-Cerenkov  part  of  the  radiation  is  attributed  to 
transition  radiation. 

In  the  present  paper,  the  transition  radiation  associated 
with  the  gas  cell  boundaries  is  included,  and  radiation  appears 
outside  the  Cerenkov  cone. 


If  the  electron  velocity  were  lower  so  that  v/c  were  close  to  but 
less  than  unity,  the  peak  in  I  would  be  pushed  to  the  left  in  Fig. 

1,  such  that  cos  0  c  =  v/c  would  be  larger  that  1.  But  the  tails  of 
the  diffraction  function  I  would  extend  into  the  physical  range  1  <_ 
cos  9  <_  -  1,  and  this  would  be  called  transition  radiation  and  be 
ascribed  to  the  passage  of  the  electrons  through  the  boundaries  of 
the  gas  cell.  Now  return  to  the  case  v/c  '•  1,  vith  the  situation  as 
shown  in  Fig.  1.  The  radiation  is  then  a  combination  of  Cerenkov  and 
transition  radiation.  The  formalism  of  Reference  12  does  admit  a 
decomposition  into  the  two  types  of  radiation,  out  is  inherently  much 
more  cumbersome. 

As  a  final  remark,  one  might  extend  the  analysis  further  in  the 
region  near  Consider  electron  bunches  emitted  from  a 

travelling  wave  Linac,  which  could  be  1  cm  long  spaced  10  cm  apart. 
Let  these  bunches  enter  the  wiggler  magnet  of  a  free  electron  laser 
(FEL) .  Then,  if  gain  occurs,  the  1  cm  bunches  would  be  subdivided 
into  bunches  of  a  finer  scale,  with  the  spatial  scale  appropriate  to 
the  output  wavelength  of  the  FEL.^3  if  the  (partially)  bunched 
beam  from  the  FEL  were  passed  into  a  gas  Cerenkov  cell,  then  the 
observed  radiation  should  be  reinforced  because  of  partial  coherence, 
at  the  FEL  bunch  frequency  and  harmonics.  This  would  lead  to  bumps 
in  the  spectrum  in  the  region  near  co  ^ . 
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Figure  1.  Qualitative  Behavior  of  the  Function  Both 

the  function  G,  from  Eq.  25  in  the  text,  and  the 
emission  angle  are  displayed  as  independent 
variables.  G1  and  G"  are  upper  and  lower  limits. 
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APPENDIX  A 


DERIVATION  OF  CERENKOV  RADIATION  FOR  A  SINGLE  FOLSE  OF 

CHARGE. 

Let  the  pulse  be  described  by 


p‘ (r, t)  *  p  q( r-  vt) 


AL 


Both  kz  and  o>  are  continuous  variables  in  this  case;  "v  is  '.gain 
along  the  z  axis.  If  we  expand  in  terms  of  a  four  dimensional 
fourier  integral, 

p'(r,t)  «  l/(2^)4  ei(wt  **  k.r)  p  '  (k,  a)  d3kd^J  A2 

It  may  be  shown  that  the  condition  Al  gives: 


p  '  (k,  to)  ■  2 tt5  ( cu  -  kzv)  p'0(k) 


A3 


where  p^(k)  is  the  three  dimensional  spatial  tranform  of  p’ 
evaluated  at  t  *  o.  All  the  fields  have  fourier  integral  rather 
than  fouries  series  expansions  and  the  energy  radiated  per  unit 
solid  angle  become 

CO  00 

r2  £  at  n  •  S  .  5  j^-2  |  £ 

A4 

l///d3r'/dt’  x  5(f't')|2 


■J 


W  (o>  ,n)  dco 


Al 


,  *  i 


The  integrand  is  a  symmetric  function  of  u  so  that. 


B<“'n)  -  I?,*  $*2M£3r'at' 


.iMtf-n-r'/e)  x 


1  u  ,  -*-.2  2 

—  ,  —  a)  vn  x  v)  M 

loir1  c 


where 


M  *  JJJJ  d3r 'dt"el(6Jt  t-n‘r  /c)p'(r,,t') 

Now  we  may  write  p'tr'^t.')  in  a  fourier  integral  representation 
p  '  (r  ' ,  t r)  =  — ^  JJJJ" d3k 'do) '  (k 1  ,o) ' )  e"*1  (a> ' * r '  ^ 


Inserting  eq.  A3  into  eq.  A7  and  the  result  into  eq  A6,  the 
integral  over  d^x'  involves  only  exponentials  and  yields 
(2ir)3  (x-  _  con/c)  ,  so  that  eq  A6  becomes 


M  *  d3lc,du,«l(a“w,')  fc,53( k'  -  fto/c) 


5  (w  ' -k  ' v)  p o  (k') 

Z  v 


Now  the  integral  over  « '  may  be  done?  because  of  the  5  function, 
o>'  is  evaluated  at  kz'v. 

M  -  J  dt'JJJd3k'eL(u~kz  v)t  6 3  (5c  *  -  a>n/c)p’(k') 


Now  do  the  integrals  over  kx',  ky‘  and  kz',  noting  that 

k  '  appears  in  the  exponential,  but  k„‘  and  kv'  do  not, 
z  * 

„  f iojt'  -iiut'n  v/c  ,  ,  ,  , 

M  =  dt  e  e  z  pa  (am  /c,am  /c,um  /c) 

J  -  x  y  z 

This  may  be  written  as 
„  f  iait'H  '  ,  a 

M  *  /  dt  e  Pa  (um/c)  a 

where 

H  ■  1-nv/c 

2  A 

If  we  let  the  time  interval  be  finite,  from  -T  to  +  T,  the 
integral  is  easily  done: 

2  * 

M  =  —  sin  oiHTpa  (noi/c)  A 

(0  n 

M2  *  4T2  sin2a)HT]  pi  (nco/c)  1  2  a 


This  result,  eq.  All  may  be  inserted  into  A5  for  cd  .  The 
factor  n  x  v  is  just  sin  9  where  0  is  the  angle  between  the 
radiation  and  the  beam  axis. 


W(o),n)  =  — — ,  -r-  m2  sin294T2  sin2a)HT |  o  '  (iW/c) 

T5S7  * 


W  is  the  energy  radiated  per  unit  solid  angle  per  unit, 
angular  frequency,  cj  .  To  proceed  to  the  total  energy,  multiply 
by  d  J2  (solid  angle)  and  integrate.  But  nz  =  cos  9  so  that,  dfl 
may  be  related  t.o  <H  : 


d  fl  =  d  (cos  0)d!2  =  -  —  dHdft 


The  functions  in  eq.  A13  do  not  contain  $  so  that,  integration 
over  ()  yields  2ir.  Thus: 


W(oi,n)dn  =  — - - a)2T2/sin29|p  o 

„  2  V  J 


2  sin  iiiHT 
(uHT)  2 


The  sin^  ojHT / (wHT)  2  factor  in  the  integral  is  peaked  at  H  =  0, 
which  by  eq.  A9  is  at  nz  *  cos  9  «*  — ,  or  the  usual  Cerenkov 
angle,  0  c.  This  function  is  more  strongly  peaked  about.  H  *  0 
for  large  values  of  T,  and  in  fact,  for  large  T  we  may  evaluate 
sin20  and  p  'Q  at  the  point  corresponding  to  H  =0.  Then  the 


integral 


2  2 

dx  sin  (ax) /(ax)  =  it/a 


may  be  used  to  evaluate  eq.  A15,  yielding 
ffwdn  *  2T  sin20c| p a  (noj/c)|  2 

A16 

The  emmlssion  was  assumed  to  occur  in  a  time  interval  from  -T  to 
+T ;  accordingly  dividing  by  2T  yields  a  rate  of  emmission,  and 
multiplying  by  v  converts  to  emmission  per  unit  path  length.  Thus 
we  obtain,  for  the  large  T  limit: 

2 

ifxdST  =  4?  ^  sin‘i9c|pe' (no>/c)l  2  A17 

where  d2E/c[ucjx  is  the  energy  emitted  per  unit  path  length 
per  unit  angular  frequency  range  cu  . 

The  corresponding  expression  for  T  not  large  is 

H" 

dw  »  ^  cudai  J  sin29|p  '  (nu/c)  |2  sin2coHT  A18 

H '  °  (ojHT)  2 

where  H"  and  H'  are  the  value  of  H  corresponding  to  9=0  and 
6  «  tt  respectively. 

Equations  Al 7  and  A18  then  describe  the  energy  radiated  per 
unit  path  length  and  per  unit  angular  frequency  range.  For  the 
non  periodic  (single)  pulse  the  radiation  has  a  continuous 
frequency  spectrum.  For  a  point  charge  q,p‘  00  is  identically  q 
and  the  usual  Cerenkov  formula  is  obtained.  Equation  A17  is 
quoted  by  Jelly,  but  only  with  the  form  factor  corresponding  to  a 
uniform  line  charge  of  length  L1 . 


APPENDIX  B 


DERIVATION  OF  EQUATION  7 

— * 

Equation  7  is  derived  for  the  case  in  which  J(r,t)  is 
expanded  in  fourier  series.  Let  the  fourier  coefficient  for  A 

be  given  by: 

T 

»(?.«)  -  7  J  at  A(?,t)  eiot  E1 

~  T  o 

Assume  that  the  green's  function  solution  for  A(r,t)  is  given: 


A(r , t) 


=  ufjjfd^'fd  f 


J(r',t) 


D(r  -  r ’ ,  t  -  t) 


B2 


where 


D(r,t)  =  4ir  5{t  "  r/c) 


Let  the  current  density  be  expanded  in  a  fourier  series: 


33 


J(r',t')  »Y  e”la) '  '  J  (r '  ,oj) 


B4 


0) 


Then  insert  B2,  B3  and  B4  into  B1  to  obtain 


t  <*,«>  >  f  |  at  eiut  fffd’r'f  df  Jj  irTTTTi  85 


«( 


t  -  f  -  |f  -  r'J/c)  ^  e”ia> ' t '  J  (r '  , 

0)' 


0)  '  ) 


B1 


dt.  ,  note  that.  t.  appears  in  the  <$  function  and  in 


Do  J 


e  ^  .  The  result,  is  t. ’  is  evaluated  at  t '  -  t  - 


A(r,<o) 


•  ~U  • ••■*///*’•'  ■*-  it?', 


Z-ioi  't  iaj  '  I  r-r  ' 
e  e  1 


/ c  -t,  i  v 

J  ,03  ) 


a) 


Do  the  integral  on  t.,  note  that 

T 


if  .  .  i  (co-co  ')!._» 

7  J  dt  «  -  6  0); 


03 


Then  do  the  sum  on  o» ' 

3 


r  i/c. 


$<*■•» jldV  ^ !<?•...  ei“i;-r'l/c 


B6 


B7 
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This  proves  the  desired  result,  B8  is  equation  7  as  used  in 


the  main  text 


Appendix  C 


TEMPORAL  STRUCTURE  OF  THE  ELECTRON  PULSE  FROM  A  TRAVELLING 
WAVE  ACCELERATOR. 

Assume  that  the  energy  of  a  single  electron  emerging  from  a 
linac  with  phase  ip  relative  to  the  travelling  wave  field  is 


E  =  Eq  cosip 


Cl 


This  relation  is  shown  on  fig.  3,  along  with  some  dots 
representing  electrons  near  the  maximum  energy  Ec,  with  phases 
clustered  about  ip  -  o  and  t b  =  2^.  Two  bunches,  separated  by  a 
phase  difference  of  2^,  are  separated  by  a  time  T-^  =  1/f Q  where 
fQ  is  the  accelerator  frequency,  which  is  fQ  =  2.85  x  10 9  Hz 
for  a  typical  S-band  accelerator  of  the  Stanford  type. 

If  a  deflection  system  with  energy  resolution  slit  passes 
only  energies  E  from  E0  to  EQ-  E  the  corresponding  range  of  phase 
Aip  is 

AE  =  E  -  E  =  E  (1-  cosAiJ; )  C2 

o  o 


For  Ai{>  small,  this  reduces  to 


AE  _  (Aip)2 

Eo  "  2 


C  1 


The  temporal  pulse  length  T2  is 


T-  =  2A  ty  T,  /2 

2  x 


C4 


or 

T2  =  Tx  •  2L']>/ 2ir 


If  C3  is  used  to  evaluate 
energy  resolution  AE/Eq 


2AE  1/2 

1  E  1 
o 


, 


1 

TT 


in  terms  of  the  fractional 


C5 


For  1%  energy  resolution,  T2/Tj_  is  about.  1/20.  The 
electrons  thus  emerge  in  short  bunches,  and  the  charge  and 
current,  when  expressed  in  a  fourier  expansion,  should  have  very 
strong  harmonic  content  up  to  and  above  the  20t.h  harmonic. 


Acoendix  D 


FORM  FACTORS 

This  section  provides  details  and  examples  of  fan  factors 
for  various  charge  distributions.  From  the  main  text,  F  differs 
only  from  P ,  the  fourier  transform  of  P ,  by  the  total  charge  q  of 
the  bunch,  so  that  for  k  =  o,  F  reduces  to  unity.  Thus  we  define 


F(k)  =  £ 

q 


p(r)eik-r 


Dl 


For  spherically  symmetric  charge  distributions,  let  k*r  =  kr 

■  ^  ^ 

where  u  is  the  cosine  of  the  angle  between  k  and  r.  In  spherical 
coordinates,  d2r  =  dudyr2dr.  Then  we  find, 


F(k) 


r 


p  (r)  sin  kr 


D2 


For  k  very  small,  sin  x  may  be  replaced  by  x  -  x  3/6  and  we 

have 


•o 

F(k)  = 

Then  the  two  terms  in  the  square  bracket  lead  to  separate 
integrals,  the  first  term  being  unity  and  the  second  is  similar  t 
the  integral  used  to  calculate  the  mean  square  radius,  <r2>, 
except  for  a  factor  k2/6.  Thus  we  have 


£  dr  r  p(r)[kr  -  k3r3/6] 

q  k  J 
^  o 


F (k)  =  1  -  k2<r2>/6 


D4 


For  a  uniform  spherical  charge  distribution  of  radius  R,  as 
well  a  a  spherical  shell  of  radius  k,  the  integral  D2  may  be 
performed  easily 


F(k)  =  3  (sin  kR  -  kR  cos  kR)  D5 

(kR)  J 

(Solid  sphere) 


F(k)  =  - 3-  sin(kR)  (Spherical  shell) 
kR 


For  a  line  charge  concentrated  on  the  z  axis,  we  may  return 
to  D1  and  let  p(r)  =  S(x)  5 (y )  p"(z),  so  that 


F ( k)  =  £  j  dz  p"(z)  e^kz  (line  charge) 
Q  J 


F(k)  =  2  sin  ( *2.) 

kZ  2 


(uniform  line  charge 
of  length  Z) 


Distorted  spherical  symmetry  may  be  said  to  occur  if  the 
scale  transformation  z‘  =  pz  serves  to  make  /o  spherically 
symmetric  in  the  prime  system.  Let  Fs  be  the  form  factor 
calculated  by  D2  in  the  prime  frame.  It  is  simple  to  show  that 


F(kx,  ky ,  k2)  =  Fs  (k^,  ky ,  k^/P) 
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